Higher order effective low-energy theories 
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Three well-known perturbative approaches to deriving low-energy effective theories, the degener- 
ate Brillouin-Wigner perturbation theory (projection method), the canonical transformation, and 
the resolvent methods are compared. We use the Hubbard model as an example to show how, to 
fourth order in hopping t, all methods lead to the same effective theory, namely the t-J model with 
ring exchange and various correlated hoppings. We emphasize subtle technical difficulties that make 
such a derivation less trivial to carry out for orders higher than second. We also show that in higher 
orders, different approaches can lead to seemingly different forms for the low-energy Ifamiltonian. 
All of these forms are equivalent since they are connected by an additional unitary transformation 
whose generator is given explicitly. The importance of transforming the operators is emphasized and 
the equivalence of their transformed structure within the different approaches is also demonstrated. 
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I. INTRODUCTION 

There is a significant recent interest in higher-order 
corrections to effective low-energy theories for a broad 
range of strongly-correlated electronic problemsiSi^ 
For example, the effective low-energy Hamiltonians of the 
Hubbard-like models contain the so-called ring-exchange 
termsi^ These terms can alter the basic properties of 
excitations;^ or shift the balance towards a new ground 
stated In addition, an accurate description of the ex- 
perimentally observed^ spectral weight transfer over the 
Mott scale also necessitates high order corrections in the 
hopping energy.-'^ Although several methods exist to de- 
rive higher-order low-energy theories, the unicity of the 
low-energy effective theory may not be obvious. In fact, 
technical subtleties appear, even in a straightforward ap- 
plication of Brillouin-Wigner perturbation theory, that 
may lead to ambiguous resultsiS beyond second order in 
the hopping. Further, unlike traditional applications of 
degenerate perturbation theory in which the eigenstates 
of the projected Hamiltonian also diagonalize the pertur- 
bation, such is not the case hereiii The unique difficulty 
that arises with Hubbard-like models is that part of the 
kinetic energy perturbation leaves the number of doubly 
occupied sites unchanged and hence must still be itera- 
tively diagonalized in the low-energy subspace. 

Low-energy effective theories play a crucial role in es- 
sentially all fields of physics. When there are large energy 
scales that are well separated from the low-energy sec- 
tor, low-energy effective theories offer the enormous ad- 
vantage of being formulated in an exponentially smaller 
Hilbert space. The price for this simplification is that 
both the Hamiltonian and the operators take a form that 
is more complicated than that of the original theory. In 
general, this form involves n-body interactions and op- 
erators that have a finite extent in space, even if the 



original theory was local. For example, making use of 
the large energy of virtual electron-positron pairs, Foldy 
and Wouthuysen derived non-relativistic quantum me- 
chanics from the Dirac equation.-^ The local gauge cou- 
pling to the electromagnetic potential in Dirac theory 
yielded three interactions, Zeeman, spin-orbit, and Dar- 
win terms, that are all non-local since they involve deriva- 
tives of the electromagnetic potentials. This was done 
with the help of a canonical transformationi^ but the 
Brillouin-Wigner type methods lead to the same result 
In the case of the Hubbard model, when the interaction 
U is large compared with hopping t, the second-order 
effective Hamiltonian in the singly-occupied subspace is 
the t-J model (including correlated hopping) where the 
spin-spin interaction is non-local in contrast to the orig- 
inal local Hubbard interaction. In higher order, n-body 
interactions, such as ring-exchange terms, also appear^ 
The appearance of n-body interactions and of non-local 
terms in effective low-energy theories is familiar in the 
Wilsonian renormalization group context^ 

There are several ways to obtain an effective low- 
energy Hamiltoniaufi^ two of which are widely used in 
condensed matter physics. The first one is the canonical 
transformation (CT) method, mentioned above, ■'■^ which 
is based on original ideas of Van Vleckii& The derivation 
of the t — J model with ring exchange^ from the Hub- 
bard model and the derivation of the Kondo model using 
the Schrieffer- Wolff transformationi? from the Anderson 
model, are well-known examples where the CT method 
was applied. A modification of that method, the con- 
tinuous CT^*-^^ and its predecessor,^*' has also found re- 
cent applications in the derivation of fiow equations. An- 
other equally popular method is the projection method, 
or degenerate Brillouin-Wigner (BW) perturbation the- 
ory. The projection method, based on ideas of KatojSi 
offers an alternative route to the derivation of, for exam- 
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pie, the t — J model with correlated hopping. ^^-^^ It has 
been also used recently in the derivation of an effective 
Hamiltonian for the pyrochloresi The resolvent method, 
based on a procedure known as Lowdin downfolding(24 is 
yet another method that is available. 

It is well known that the BW and the CT methods give 
equivalent results at low (second) order, as exemplified 
by the t—J model with correlated hopping. Although the 
equivalence of these approaches in higher order might be 
not obvious, one intuitively expects that they are. This 
has been supported by Kleini^ who has given a formal 
proof of the equivalence between various forms of degen- 
erate perturbation theory. Generally, in the presence of 
a well defined small parameter, the perturbative expan- 
sion in powers of such a parameter should not depend on 
the method employed. In other words, an effective low- 
energy theory can be presented as an n-th order power 
series in the inverse of some large energy scale {t/U m the 
Hubbard model), and all methods are expected to yield 
equivalent forms of the theory up to that same order n. 

In this paper we show in some detail how CT, resolvent 
and BW methods can be applied to the Hubbard model 
to obtain equivalent effective low-energy Hamiltonians up 
to fourth order, that is 0{t'^/U^). By equivalence, we 
mean that all three methods yield Hamiltonians that are 
related via a unitary transformation. Although, in order 
to be specific, we work with the Hubbard model, it will 
be clear that the procedure can be trivially extended to 
other models, including models that involve expansion of 
pure spin models about the Ising limit. Several new is- 
sues appear in deriving higher-order low-energy effective 
theories. First, one may find amusing that the low-energy 
Hamiltonians obtained from three different CT methods, 
one of Refs. |^|^ one of Ref. and the one introduced in 
this work, appear to be different in each case. However, 
we show that they are all connected by an additional uni- 
tary transformation that leaves the block-diagonal form 
invariant,?- That is, this unitary transformation converts 
these different Hamiltonians one into the other. Second, 
in the case of the BW method one should be careful in 
dealing with: (i) the orthonormalization of the projected 
eigenstates, since the latter are not necessarily orthogonal 
even if the initial basis is, (ii) the energy dependence of 
the expansion, since the energy should be evaluated iter- 
atively using results from previous steps, in the spirit of 
Rayleigh-Schrodinger non-degenerate perturbation the- 
ory. The final low-energy theory again appears to be dif- 
ferent from that obtained from the CT approach of Ref. 
0, but a unitary transformation within the low-energy 
subspace shows their equivalence. Third, in the so-called 
resolvent method, which is similar to the BW approach 
but works more directly with the Hamiltonian matrix 
rather than with the eigenstates, one needs to perform 
the orthonormalization of the eigenstates in an iterative 
procedure. For the derivation of the t-J model with cor- 
related hopping within this method see, e.g. Ref. 
The fourth-order low-energy theory obtained with this 
method is also unitarily equivalent to the results of the 



BW and CT approaches. 

Finally, as in Ref. |^ we emphasize that in order to 
compute correlation functions or spectral weight within 
the low-energy theory it is important to transform the 
operators corresponding to observables. The omission of 
such a transformation (see, for example, the recent work 
using the BW method for the pyrochlores^'^^) has to be 
explicitly addressed. We also demonstrate the equiva- 
lence of the transformed structure of the operators within 
the different methods. 

Our paper is organized as follows. Sec. Hll introduces 
the model and notations. Sec. IIIII describes CT methods 
and introduces the unitary transformation that allows us 
to demonstrate the equivalence of the various low-energy 
effective theories. Sec. II VI is devoted to the BW method. 
Sec. |V] describes the resolvent method and finally Sec. 
IVII presents a generalized transformation from which all 
perturbative results can be derived. In Sec. IVIII we give 
a brief discussion of relevant experiments. We conclude 
with Sec. IVIIII Appendix 1X1 contains miscellaneous re- 
sults. 



II. MODEL 

We consider the Hubbard model, conveniently written 
in the form^ 

H = To + Ti+T^i + V, (1) 

where the 0{t) kinetic energy operator in second quan- 
tized form has been divided into three terms by using 
projection operators. The first operator, Tqj includes 
the projection operators that ensure that the number of 
doubly-occupied sites does not change because of hop- 
ping. It is precisely this term that complicates the appli- 
cation of traditional degenerate perturbation theory to 
the Hubbard model. The projection operators included 
in the kinetic energy operators Ti and T-i make sure 
that these operators increase or decrease the number of 
doubly-occupied sites by 1, respectively. All of these 
terms are proportional to the hopping matrix element 
t. Note that generally speaking, there can be hopping 
matrix elements to arbitrary neighbors but we take all 
these terms to be of the same order in the expansion pa- 
rameter t. The Hubbard on-site repulsion, written as V ^ 
is proportional to U . 

These notations implicitly use the classification of the 
Hilbert space in subspaces with different numbers of 
doubly-occupied sites. Namely, every eigenfunction of 
the Hubbard Hamiltonian ^ can be split in the series of 
orthogonal pieces 

oo 

IV') = E 1^™) = 1^") + i^i) + 1^2) + . . . , (2) 

where the subscript enumerates the states with m = 
0,1,2,... doubly-occupied sites. Evidently, the operator 
V is diagonal in this basis and has the eigenvalue mU in 
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a state containing m doubly-occupied sites. For the rest 
of the paper we set U = 1, and consider t itself as a small 
parameter. When a clarification is needed we will restore 
the actual t/U dependence of the expression. 

At t = 0, the solution of the eigenstate equation is 
simply a set of highly degenerate states separated by en- 
ergy U. It is assumed that U is much larger than the 
bandwidth (of order t) so that, with t finite, states clus- 
ter around the values mil and are separated from each 
other by a Mott gap where no states occur. In other 
words, we assume that switching on t does not lead to 
a crossing of levels between the lowest energy manifold, 
mU with m = 0, and all other manifolds with m > 0. 

Evidently, many Hamiltonians can be cast in the form 
of Eq. (QJ . It suffices to have a large term in the Hamil- 
tonian that can be easily diagonalized and leads to states 
that are separated by some large energy scale U. Other 
terms in the Hamiltonian that couple these original states 
and are proportional to some small parameter can be 
denoted by TAm where Am indicates that the operator 
couples states separated by Amll. Although we consider 
only Am = — 1, 0, -|-1, it will be obvious how to generalize 
our proof to more general values of Am. 

III. CANONICAL TRANSFORMATIONS 

The CT method is probably the most commonly used 
method to find low-energy effective theories. We briefly 
recall the known results for the problem at hand and then 
move to another version of the CT approach that gives a 
result that is unitarily related to the first one. 

1. Method 1 

To derive the higher-order effective Hamiltonian for 
the Hubbard model, the CT method was applied in Refs. 
1^0. There the t^- and i^-order Hamiltonians were ob- 
tained for the half-filled case and i^-order Hamiltonian 



was found for an arbitrary doping. We will simply repeat 
the basic idea and the results obtained in Ref. i^. The 
effective Hamiltonian was obtained from the Campbell- 
Baker-Hausdorff expression 

rfj^ ^e-^He^ ^H+[H,S] + ^[[H,S],S] + ... , (3) 

where the generator of the transformation S is truncated 
as 

^ ^1 -t- ^2 + ^3 + . . . 

where 5„ oc t". The role of each Sn in this series is to 
eliminate the corresponding t"-order off-diagonal terms 
in the Hamiltonian in Eq. (jS)) , which change the number 
of doubly-occupied sites. By assumption however, S does 
not contain terms that preserve the number of doubly- 
occupied sites. The remaining freedom to perform a uni- 
tary transformation within the singly-occupied subspace 
will be discussed later in this section. 

Given the explicit form of the Hubbard Hamiltonian, 
Eq. (Q, one readily finds that 

Si - r_i - Ti . (4) 

Using S* = and keeping terms up to 0[t^) in Eq. (0) 
the t — J model with correlated hopping is obtained. To 
derive the higher-order Hamiltonians we have to truncate 
S at higher order and determine the operator expression 
of SnS required to ehminate the off-diagonal terms of 
the i" order. These off-diagonal terms are generated by 
the commutators of H with Sn' 's in the previous orders. 
Generally, one needs n— 1 terms in S to obtain the theory 
valid to the order i". We list here the generator iS'2 for 
completeness 

^2 = ToT-i - TiTo - r_iTo + ToTi , (5) 

and simply reproduce the result of the procedure de- 
scribed above carried out in Ref. 6 to the f^-order. 



-T^iT^Ti + T-iToTiTo + ToT^iToTi - ^{T^iTiT^ + T^T^^Ti) . (6) 

I 



The value of V is taken to be zero since we are in the 
singly-occupied subspace. In deriving this expression the 
identity T_i|-0o) = 0, where 1-00} is any of the singly- 
occupied states, is used. 

We would like to add an interesting technical detail 
to the discussion of this method. The t^-term in the 
generator, 5*3, although necessary to eliminate the t^- 



order off-diagonal terms in Eq. (jSjl, does not contribute 
to nlf^'^^^ explicitly. That is, the 4i/i order effective 
Hamiltonian of Eq. lO can be obtained using S = S'i-l-5'2 
only, simply neglecting the remaining off-diagonal terms. 
Similarly, the diagonal terms in the 3rd-order Hamilto- 
nian Ti-^lj are all generated by 5*1 alone. This is because 
the original Hubbard Hamiltonian does not contain any 
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"bare" off-diagonal terms of order higher than t. All such 
higher-order off-diagonal terms are the result of commu- 
tations of T-terms. This shows some additional internal 
structure of the model. Since we will need the generator 
5*3 for the discussion of the transformed operators but it 
was not written out explicitly in Ref. we present it in 
Appendix El Eq. (|XT|l . 

The expression for the Hf^jF"^^ in Eq. © seems to 
have all possible combinations of Tq, T-i, and Ti, except 
for one "missing term" : TqT-iTiTq. Although there is no 
general principle which would require presence of such a 
term in the effective Hamiltonian, its absence makes one 
curious about its whereabouts. The fate of this term will 
be clarified in Sec. IIII 31 

2. Method 2 

In the same spirit, a different way to formulate an effec- 
tive theory using the CT approach is to apply consecutive 
unitary transformations 

n'^^^ = . . . e-^^e-^^-e"^^ H e^^e^^e^^ .... 

We call this the "consecutive CTs" approach. The idea 
for each Sn is to eliminate the off-diagonal terms of the 
nth order remaining from the previous, n — \ order CT. 
In each of the transformations the expansion formula Eq. 
Q is applied and all terms up to a desired order in t 

are kept. Thus, after Ui = e^^ is applied to the origi- 
nal Hamiltonian the off-diagonal terms of order 0{t) are 
eliminated and show up only in O(t^). The next trans- 
formation moves the off-diagonal terms to 0{t^), and so 
on. Generally, the generators Sn in this approach are 
different from the ones in the previous approach, that is 
Sn ^ Sn- However, one can check that for the Hubbard 
model generators Si = Si and S2 — S2, Eqs. Q and 
©. Note that S3 is indeed different from S3. 

We would like to remark that the derivation of the 
higher-order effective Hamiltonian within this "consecu- 
tive CTs" approach is more straightforward than the ap- 
proach of Ref. 6. Also, the absence of the S3 contribution 
to '^e^-j is much less enigmatic here. Namely, since 53 is 
the off-diagonal order operator, its only commutator 
which can give 0{t^) contribution to the Hamiltonian is 
with the off-diagonal operator of the order t. However, 
after e'^^ is applied the only operator of order t remaining 
in the transformed Hamiltonian is the diagonal operator 
Tq. This removes the need to know the explicit form of 
S'3, although it is formally still necessary to eliminate the 

off-diagonal i'^-order terms in '^e^j- 

Surprisingly, the final 4:th order result of the "consec- 
utive CTs" approach is different from Ti.f^^j!'^'^^ Eq. ©: 



The difference concerns the above-mentioned "missing 
term" TqT^iTiTo and the terms T^T_iTi and T^iTiT^. 
In fact, in this version of the effective theory the original 
"missing term" is found, but the two analogous terms are 
missing. This "mystery" is unveiled below. 

3. Additional unitary transformation 

Let us first make the following observation. One can 
consider the following unitary transformation: 

K// = e-^"Hf^7e^« , (8) 

with the generator 

So = 7(ToT_iri - r_iTiTo) , (9) 

where 7 is a real number and plays the role of an "angle 
of rotation". Note that such a generator: ('ij is explicitly 
anti-hermitian, (ii) is diagonal (does not change the num- 
ber of doubly-occupied sites), (Hi) is 0(t'^), and (iv) is 
real. The operator Sq is the lowest order operator satisfy- 
ing (i)-(iv) which one can construct using Tq and T±i. It 
is also the only operator of such kind in the 0{t^) order. 
Therefore, the only contribution from such a transforma- 
tion to Ti-'^ff will be from the [Tq, Sq] commutator and it 
will generate additional terms of t^-order of the form, 

SUeff = l{T^T^iTi + T^iTiT^ - 2ToT^iTiTo). (10) 

As we will see, the multiplicity of Hamiltonians that arise 
once the high energy scale is eliminated all differ by the 
terms appearing in STief / • Choosing the "angle of rota- 
tion" 7 = 1/2 and applying the transformation Sq to the 

nf^f^\ Eq. ®, one readily obtains H^ff^^ Eq. 0. 
As a result T_iTiTq -I- H.c. are replaced by the "missing 
term" TqT^iTiTq. Clearly, different choices of the "angle 
of rotation" will give different fractions of those terms in 
the result. In fact, a recent study, Ref. '25', used a contin- 
uous CT approach to the Hubbard model and obtained 
an effective Hamiltonian which would be equivalent to 
the choice 7 = 1/4 in Eq. 

Thus, the Ath-order effective Hamiltonian for the Hub- 
bard model can take an infinite number of unitary equiv- 
alent forms, all connected by the transformation in Eqs. 
©, ©■ AH these models possess the same energy spec- 
trum and correlation functions and thus are equivalent. 
From this point of view, the reader should not be sur- 
prised when, in the next section, we find that the BW 
method gives a result that is different from Eq. 10. 

Although the unitary equivalence of the models is 
a rather natural property, it is certainly unfamiliar in 
lower-order effective theories. Furthermore, such a uni- 
tary equivalence should be common to all higher-order 
theories. As the order of the perturbation theory is in- 
creased the number of block-diagonal, anti-hermitian op- 
erators one can construct will also grow, providing one 
with a broader variety of unitarily equivalent forms of 
the effective Hamiltonian and corresponding operators. 



5 



4- Operators 

It is important to note that in a low-energy effective 
theory all operators should be transformed along with 
the Hamiltonian. Then the expectation values of the ob- 
servables can be calculated in the singly-occupied mani- 
fold. The transformation is different depending on which 
canonical transformation method is used. Using the first 
CT above, the standard expression for the transforma- 
tion is 

O ^ e-^Oe^ ^0+[0,S] + ^[[H, 0],0] + ... . 

Again, we consider as an example the operator Oi , which 
increases the number of doubly-occupied sites by one. 
Using S = Si + S2 from Eqs. Q and (jSJ and utilizing 
the property T^i\iPq) = we obtain, to the order 

Oi = -T_iOi + (T_iTo - ror_i)Oi, 

which coincides with the expression we will obtain with 
the BW method, Eq. (|34|l . To obtain the next-order ex- 
pression for the transformed operator one needs to know 
an explicit expression for the generator S3 (see Appendix 
IXI Eq. (jAljl l. Using it, some algebra reveals that, to or- 
der 

Oi = -r_iOi + (T_iTo-ror_i)Oi 

+ ^T_iTiT^iOi - ^T\TiOi (11) 
-T^T^iOi - T^iT^Oi + 2Tor_iroOi 

which should be compared with the result of the BW 
method given in Appendix 1X1 Eq. (|A2|I . 

IV. BRILLOUIN-WIGNER METHOD 

We proceed to show, up to order 0{t'^) (= Olf^/U^)) 
for the Hamiltonian, that degenerate BW perturbation 
theory can be organized in the spirit of the Rayleigh- 
Schrodinger (RS) perturbation theory to lead to the same 
low-energy effective theory for Hamiltonian and opera- 
tors as the consecutive CT method. We will take a de- 
tailed approach that shows all the subtleties. 

Let us consider one eigenstate of the full Hamilto- 
nian ^ with the eigenvalue E. It obeys the Schrodinger 



equation 

{E-To-V-Ti- T_i) \iP) = 0. (12) 

Although we have not explicitly written quantum num- 
bers for E and \4>), we have to remember that we have a 
matrix equation with many eigenvalues and correspond- 
ing eigenstates. We look for the effective theory that 
describes the states that evolve from the lowest energy 
sector, TO = 0, taking into account virtual excitations 
into TO > states perturbatively. One can also write ef- 
fective theories that are valid for any of the subspaces2& 
with TO > 0. 



1. BW expression for IV') 

We would like to rewrite the eigenstate {ip) in a way 
that will allow us to take into account higher-energy sec- 
tors with TO > 1 through an iterative procedure. Let Q 
be a projection operator that removes all components of 
that are in the to = (singly-occupied) subspace. We 
have [Q,E — Tq — V] = since Tq + V does not change 
double-occupancy. One can find then an iterative expres- 
sion for Q\ip) directly from the Schrodinger equation Eq. 

m 

Q\^)^^-^^—^Q{Ti+T^i)\ij). (13) 

Inversion of the operator E — Tq ~V does not cause any 
problem when there is a Mott gap since the denominator 
has only non-vanishing eigenvalues. Indeed, E — Tq is at 
most of order of the bandwidth (proportional to t) while 
the operator Q ensures that the smallest value that V 
takes is U. The complete eigenvector \tjj) has components 
in the to = subspace \ipo) = {1 — Q) \if}) that we need 
to determine. We assume that |?/'o) is a member of an 
orthonormal subspace (i/'gl'^/'o) = ^%.ipo- The subscript 
to a ket means that it has components only in the 
TO = subspace. This procedure leads to the standard 
BW expression for \tp)^ 

IV') = l^o) + ^_^^_^ Q(ri + r-i)|^), (14) 

which can be solved perturbatively by iteration. Using 
T-ilV'o) — and the fact that we cannot come back to 
the TO = subspace in any of the intermediate steps 
(QT_iT'i|V'o) — 0) we find, iterating Eq. ((Tljl three times. 



1 

E-To-V 



1 



E-To-V 



Ti 



1 



E-Tq-V 



Ti 



(15) 



E-n-v 



E-Tn-V 



Ti 
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We took into account the projection operators Q so that 
the above equation contains only the terms for which Q 
equals unity. The above equation Eq. (fT^ generates the 
usual Brillouin-Wigner perturbation theory. One recog- 
nizes that the second term and the last term in Eq. H15() 
are components of the eigenvector in the subspace with 
m = 1 doubly-occupied site li/ii), while the third and 
the fourth terms are components with m = 2 (|'02)) and 
m ~ i (IV's)) doubly-occupied sites, respectively. This 
form (without ji/'s)) suffices for our derivation of the ef- 
fective theory to order 0{t'^). With this effective Hamil- 
tonian, one will be able to find the component of the 



eigenstate in the singly occupied subspace |'0o)- Given 
IV'o), all the components of the eigenvector \'ip) in the 
m = I and m = 2 subspaces are already completely de- 
termined by Eq. (|15|l . 

The subsequent treatment of Eq. (|15|l to generate a 
low-energy theory is the following. The denominators 
{E — Tq — V)~^ are not singular because they correspond 
to the energy in the bands with m > 0. Using V S> 
E ^ Tq, one has to expand the energy denominators in 
Eq. (|15|l to the required order in t. Let us list here the 
results of such an expansion of Eq. (|15|l order by order 
{U = 1). We find 



1^(2)) 



= (1-Ti)|v4'^ 



= [l-n-{E-To)T, 



1 - Ti + TqTi - TiE - -T^iTI - TiE^ - T^Ti + 2TqTiE 



(16) 
(17) 

(18) 
(19) 



where the superscript describes the order of approxima- 
tion, that is |'0o ) 1^ the component of the third-order 
eigenstate |-)/)(^)) in the singly-occupied, m — subspace: 
1^/;^^)) = (1 - (3)|i/'^^^)- Note that, as usual, the order 
j-n-i -j-j^g expansion of the eigenvector correspond 
to the order in the matrix elements of the Hamiltonian 
and of the corresponding energyi^ That is, one computes 
the first-order i?^^) using zeroth-order eigenstates 
the second-order theory _ff is formulated with the first- 
order basis I'i/''-^-'}: etc. Thus, for the fourth-order effective 
Hamiltonian we will need Itp^^^)- Another detail concerns 
the explicit dependence of IV'''^'') and in Eqs. (|18|l . 

(|19|l on the energy E. This issue will be resolved later by 
evaluating E in iterative manner. 

One can see that in all orders the full eigenstate |V') 
is built from the m = states Itpo) by including the 
off-diagonal transitions to the upper-band. The zeroth- 
order approximation Eq. H16|l corresponds to tak- 
ing only the first term in Eq. (|15|l and neglecting all 
the upper-band excitations. For the first-order approxi- 
mation we needed to include the second term in 
Eq. ((TB|l . Using that E — Tq will be of the order of t we 
have simply neglected E — Tq in the denominator which 
yielded the first-order state, Eq. ifTTI) . The second-order 
state IV'''^'') is obtained from the first three terms in the 
BW series Eq. (I15II . In the third term we can neglect 
E — Tq again, but the second term needs to be expanded 
in E — Tq. This leads to the second-order state given in 
Eq. (|18|l . A term of the form /2, which does not con- 



tribute to the effective Hamiltonian to that order, was 
dropped although it can appear in certain observables. 
To obtain the effective Hamiltonian valid to order t'* one 
has to generate the third-order wave function j^^^^) from 
Eq. (|15|l . This is obtained from the contributions of all 
terms in Eq. (|15|l , although the fourth term cx Ti can be 
neglected since it only contributes to the effective theory 
to order t^. The expansion of the denominators provides 
us with the E-, and i?^-dependent terms in the result- 
ing expression for IV'*-'^^), Eq- (|19|l . Let us emphasize 
here that in the _E-dependent terms of the second- and 
third-order eigenfunctions, one should expand E as well, 
and keep only the terms to the required order. Thus, for 
\ip'''^^} we will only need the expression for E that is valid 
up to order i, while for order is required. 

We also point out that equations relate the 

"full" eigenstate j?/') to the the state in the "projected", 
m — subspace. Since our goal is to have a low- 
energy theory which operates with the projected \ip()) 
states only, one should take into account the fact that 
although the eigenstates are orthonormal ('0'IV'') = S^',ip 
this is not true for the m = component alone, that is 
(tAoI^Ao) 7^ ^Vo,'/'o- "^^^s is evident from Eq. Q. There- 
fore, if one wishes to obtain an effective theory that takes 
the usual Hermitian form, one needs to orthonormalize 
the projected basis to the required order in t. This dif- 
ficulty will appear in higher order, but let us first repro- 
duce the well-known second-order results. 
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2. Derivation of the t-J model 

We begin with the zeroth-order approximation |-)/)(*')) 
given in Eq. (|15|l . Since the eigenfunction involves only 
the m = subspace, the diagonal part of the Hamiltonian 
(0 solely contributes to the eigenvalue equation 

= (^rirol4°^) = ES^,^^^ . (20) 

Therefore, Ti-l^jj: = Tq is our effective Hamiltonian to 
order t. We rewrite Eq. (|20() . for future reference, as 

Wl,"^) = E\4'^) . (21) 

To next order, we need the first-order state IV'^^'') from 
Eq. (fTTji . One can see that the states \ipQ^^) in Eq. l(T7| 
are still orthonormal to order t, that is [^''^^^iIj'-^'>} = 

(tAq IV'o^^) + 0{t'^) = (5^',^. With these states one ob- 
tains 

(^''(DliJIV'W) = (V-^^i^lTo-T^iTilV^W) (22) 

where we used that E ~ 0{t) and that with the required 
accuracy we can neglect 0{t^) term in the right-hand side 

and use orthonormality to replace (V'o IV'o^^) by S^'^^^g. 
It is clear that the content of the brackets in Eq. 1)22(1 
plays the role of an effective second-order Hamiltonian 

nff\ = To - T^iTi . (23) 

Writing this result in terms of second-quantized operators 
and recalling that J = At'^/U in the second term, one 
recovers the t — J Hamiltonian with correlated hopping. 
Again, for future reference, we write: 

(To-T_iri)l4")=^^l4'^) ■ (24) 

Contrary to the non-degenerate perturbation theory, 
the states in the m = subspace change as we improve 
the approximation. There is a link between the states at 
various order as will be discussed in more detail in Sec. 
ITvel below. 

3. Transformation of operators in the t-J model 

Let us pause momentarily to develop transformation 
rules for the operators that should be used at this level 
of approximation. Consider, as an example, an operator 
0_i that decreases the number of doubly-occupied states 
by one. Naively, one would expect that it has zero expec- 
tation in the case of the t-J model that is defined in the 
singly-occupied subspace. The correct way to proceed is 
to notice that the matrix elements of 0_i in the basis of 
first-order eigenstates Eq. (|17|l are given by 

(^'W|0_i|V'(^)) = (V^ri - 0_iTi|4^)) . (25) 



Since the above expression is valid for any eigenstate, to 
first order the effective operator 

O^ll = -0-iTi (26) 

should be used to compute any matrix element of the 
original operator solely in terms of projected eigenstates 
\tjjo). Both the Hamiltonian Eq. (|23|l and the operator 
Eq. (1211) coincide with the CT result. 



4- Third order 

There are new technical issues that appear in orders 
higher than second. The first problem is that since 
(V;'(2)|V;(2)) = )|1 + r_iri|4'^) + 0(i3), the second- 

(2) 

order eigenvectors in the singly occupied subspace IV'o ) 
in Eq. H18|l do not form an orthonormal set to order 
t^. This issue is easily resolved. To obtain an eigenvalue 
problem in standard form we define the orthonormal ba- 
sis \ip'^^) by 

l4')) = (l + T_iTi)-'/'|^[,'^) (27) 

where the square root needs to be expanded to order 
to give 

IV^')) =(^l-T,-{E- To)T^ - ^T^iTi) (28) 

This second-order eigenvector IV''-^-') should be used to 
obtain the third-order effective Hamiltonian. This is 
where one encounters the second difficulty. It concerns 
i?-dependent terms that have to be treated carefully. 

Suppose IV'*-^-') in Eq- H28|l is an eigenstate with en- 
ergy E. Since the order of the £'-dependent term should 
match the i^-order of the eigenstate IV'*-^''), E should be 
expanded in t. Clearly, we can replace E appearing in 
Eq. H28|l by its first-order term in powers of t. To this 
end we write 

i^l^[,'^)=^3>o") • (29) 

(3) 

Since we can anticipate Hg^j to have the form 

n^^^f = To - T_iTi + 0{t^) , (30) 

we can use E\ip^q'^) = To|(^[,^^} + 0{t^). This is an it- 
erative procedure and it can always be justified by the 
consistency of the result with the initial expectation out- 
lined in Eq. The replacement E\(p^^'') = To\ip''Q^) 

(2.) 

is valid for any of the eigenvectors \(pQ ) of the effective 
Hamiltonian. After making this change in Eq. H28|l . the 
second-order eigenstate reads 
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IV'^'^) - (^1 - Ti - TiTo + ToTi - if-iTij Wf^). (31) 

This procedure allows us to find a Hamiltonian matrix. We will show in Sec. II V 61 that there is a relation between 
eigenvectors and eigenvalues at different orders. With the above state l^/''-^-'}, the eigenvalue problem takes the form 

(V^'^^^liJl^^^)} = - T^iTi - iT^iTiTo +T_iToTi - ItoT^iTiI^')) = , (32) 

which is clearly consistent with our expectation for the eigenvalue E expressed in Eq. H3U|) . As in the lower orders, 
the effective low-energy Hamiltonian can be directly read off this equation. Thus, up to order 0{t^), the effective 
Hamiltonian is given by 

nfjf = To - T^iTi - ^T^iT^To + T^iTqTi - ^ToT^iTi , (33) 

which agrees with the CT approach to this order (see Sec. Ifff|l . 

The last equation, when taken out of context, may suggest that it contains terms that are not allowed within BW 
perturbation theory. Indeed, the general expression with projected wave function Eq. (|14(l shows that one cannot 
come back to the singly occupied state in any intermediate state. The series should contain only "proper" terms. 
So, terms like T_iTiTo in the above expression leave the impression that they should be forbidden since Tq does not 
change the double occupancy. However, our derivation clearly shows that these terms come from expanding E using 
lower order results. Hence they are in fact perfectly acceptable. 

The matrix elements of any operator to second order should be computed with the state l'^'-^-') from Eq. H31|l . For 
the operator Oi, which increases the number of doubly-occupied sites by one, we obtain 

(V'(')|Oi|V'(')) = (^f 'I -T_iOi + (T_iTo -ToT^OOil^"') , (34) 
that also coincides with the CT result for the transformed operator, as discussed in Sec. IHII 



5. Fourth order 

To obtain the effective Hamiltonian valid to the order one needs again (i) to orthonormalize the states IV'o'^''), 
Eq. H19I) and (ii) to transform the £'-dependent terms to equivalent operator expressions found in the previous steps. 
The procedure outlined for the second-order should be followed here. From the normalization condition of 

in Eq. ()19fl . one finds that iV'o'^'') should be orthonormalized with the help of 

IV-^'^) = (1 +T_iTi - 2r_iToTi + 2T^iTiE)-'/^ \ipl^^) . (35) 

The subsequent expansion of the square root to order is needed. The resulting expression for the eigenstate 

will contain a term —TiE. One has to to replace E\(pQ ) with (Tq — r_iTi)|iy9g }, which gives E to second order in t 

when liySg^'} is an eigenstate. This leads to the second- and third-order termSi^ In all other i?-dependent terms in Eqs. 

((T^ . l|H5|) we can use E\ip'^'') ~ To\(p'^^) as before since this already produces the terms of order t^. Straightforward 
algebra finally leads to 

IV'(^)) = (^1 - Ti + ToTi - TiTo + - lr_iTi - ^T^iTf - TiT^ - T^Ti + 2ToTiTo (36) 
-ItiToTi + ^T^To - ^ToTi^ + T^iToTj - T^jTiTq + ^TjT^iTi^ |<^[,')), 

where |(^o^'') is the orthonormal set of the m — 0, singly-occupied eigenstates. With the help of this form of IV'*''^'') the 
third-order expressions for the matrix element of an operato r (Oi) can be obtained (see Appendix 1X1 Ea. lfX2|l). 



Rather cumbersome, but still straightforward calcu- 



lations for the eigenvalue problem with I'lp'^^^) yield the 
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fourth-order effective Hamiltonian 

^eff ~ ^eff 



(37) 



that is identical to the one obtained by the consecutive 
CT method, Eqs.Q. This effective Hamihonian, within 
a unitary transformation in the singly occupied subspace, 
is the t-J model with ring exchange and various corre- 
lated hoppingsi^ 



6. Connection between eigenstates at different orders 

There is a connection between eigenvectors in the m — 
subspace at different orders in t. For definiteness, let 
us consider the eigenvalue problem defined by 'H^^fj: in 

Eq. ijSSl)- If Iv'o^) is an eigenstate of Ti-fff, then the 
corresponding energy to order t'^ is given by 

(To - T^iTi + Oit^)) |<^(')) = E\cp^^'>) . (38) 

We wish to rewrite this equation as 

(To - r_iTi) |^(^)) = {E- Oit')) Ivp*'^ (39) 

Note that the left-hand side of this expression contains 

(2) 

Tq — T^iTi = 'H^jp the effective Hamiltonian at the 



(2) 

previous order. The eigenvalue problem for n^j/j: gives 

eigenstates at the previous order, \(p'^^) (= |'0o^^) in Eq- 
(|f 7|l ). Taking the hermitian product of Eq. (|39|l with an 

eigenstate ((ySg | we have, by applying Tq — T^iTi on the 
bra, 



(2)\ _ 



V(l)|,o(2)\ 



(40) 



From this, one concludes that to order i^, 



(i?-£;')(4'V^'^)-o 



(41) 



Hence, to that order, either (vo ^"^ I'/'o^^) = or E' = E 
or both. The generalization of this result means that an 
eigenvector in the low-energy subspace can have a non- 
zero overlap to order t" with an eigenvector of the order 
i"^^ theory if and only if the energies agree to order 

(n) 

Note, by the way, that the diagonalization of Hg^j will 
in general give us energies that contain all powers of t. 
Nevertheless, E will be valid only to order since the 
higher orders can be modified when the matrix elements 
of Hef f are calculated to higher order. A degeneracy that 
exists at a given order in t can be lifted at the next or- 
der. Energy levels of different symmetry can cross when 
evaluated at different orders in t so that the ground state 



of, for example, the t — J model with correlated hopping 
can be different from that of the model that includes ring 
exchange. 

V. RESOLVENT METHOD 

Another approach to deriving the low-energy effec- 
tive theory is the resolvent method, ft is based on 
an iterative execution of a procedure known as Lowdin 
downfoldingi^ It bears a lot of similarity with BW per- 
turbation theory but works more directly with the Hamil- 
tonian matrix rather than with the eigenstates. Of all the 
approaches considered in this paper the resolvent method 
requires the least amount of algebra. 

We start by defining the projection operators Pm, m — 
0, 1, 2 and P>. Operator Pm projects on a subspace with 
m doubly occupied sites. Operator P> projects on a 
subspace with more than two doubly-occupied sites 



(42) 



m=0 



For our purposes it is convenient to rewrite Eq. Q as 



(43) 



where i,j = 0,1,2, and ">". The eigenvector equation 
can then be written in the following block form 



To T_i 
Til + To T_i 

Ti 2 + ToT_i 




E 



( IV^o) 

IV^i) 
Vl^>), 



(44) 

where U> ^ P^HPy. Similarly to Eq. lO in the BW 
formalism one needs to keep only terms up to \^p2) to 
derive the effective theory to order t*. Thus, Eq. 144(1 
should suffice for our goals. It is convenient to rewrite 
this equation by components 

E\ijo) ^Wo)+T,M, (45) 

E\^i) = (l + To)|V^i)+Ti|V^o)+T_i|V^2), (46) 

E\iP2) = (2 + To)|^2)+Ti|V^i)+T_i|^>), (47) 

£;|V>) =^>|^>)+Ti|^2). (48) 

Now we eliminate all components of lip) one by one, start- 
ing with !?/;>) until only lipo) is left. From Eq. (|48l) we 
obtain 

|V'>)-(i?-W>)"'Ti|^2)=0(t)|V2), 

where we take into account the fact that the operator in 
brackets is non-singular and is 0(1). This expression for 
V'>) is substituted in the equation (|17|l for |?/'2) to give 
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1^2) ^{E-2-To- 0{t)) VilV^i) = (^-^Ti + 0(i2)^ l^i), 



(49) 



where we expanded the denominator and kept terms of order t since the higher order terms do not contribute to the 
theory of the required t'^ order. This latter equation is used to ehminate [^2} from Eq. (|46|) . Thus we have 

1^1) = [-Ti -{E- To)T, - \t_,TI -{E- nfT, + 0{t^)^ |^o), (50) 

where again the expansion of a denominator has been performed to the required order. FinaUy, we obtain an equation 
for IV-o) 

E\^o) = (to-T_,T,-T_,ET,+T_,ToT,-^T\T^ 



- T^iE'Ti - T^iT^Ti + 2ET^iToTi + 0(i^)j IV'o)- (51) 

This is not the "true" eigenvalue equation since it contains _B-dependent terms in the right-hand side, similar to the 
BW case. We rewrite it then by transferring £'-dependcnt terms to the left as 



(52) 



where we omitted i^-order terms. Below we will assume that operator expression accuracy is up to O(i^). 

The derived equation for 1-00) still does not have the form of the Schrodinger equation E\il}) = TL\i)). To transform 
(|52|l into a Schrodinger equation an orthogonalization procedure similar to the one we used in BW calculations must 
be performed. We introduce |xo) such as 



Substituting in Eq. I|52|l . left multiplying by a and then expanding the square root keeping third order terms, the 
left-hand side is E\xo) and the right-hand side still has some iJ-dependent terms. Transferring them to the left we 
find that |xo} satisfies 

e[i + \ {T-iTiTo + ToT^iTi)^ |xo) = - T^iTi - ^T^iT^To - ^ToT^iT^ + T^iT^T^ (53) 

+T^inTin+ToT^iToTi + {T^iTi f - ^T\T^ - T^iT^T^ |xo). 

The left-hand side of this equation still does not have the desired form. An extra orthogonality transformation 
analogous to that performed with a above is required 

1 1 ^ 



Ixo) =\l + ^T^iTiTo + -ToT^iTi j |(po), 

with the subsequent square root expansion. The resulting eigenvalue equation for j^^o) by the resolvent method finally 
takes the Schrodinger equation form E\ipo) — Ti^fj^^'lipo) with the effective Hamiltonian given, to fourth order, by 

Hfff = T^yf -I- ^{T^iTiT^ + T^T^iTi) - ^ToT^iTiTa . (54) 



r 



This effective Hamiltonian does not coincide with the Eqs. ©,07 and (|37|) . It is however unitarily related to 



other forms of Tigt''* we have obtained so far, namely all others through the transformation Eqs. @, 110|) . For 
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example, the "angle of rotation" 7 = 1/4 transforms Eq. 

(EH) back to n^^ff^^ in Eq. ®. 

We note that transformation of operators can also be 
devised within the resolvent approach in a manner sim- 
ilar to the BW calculations, Sec. IIVI using the above 
relations between I^Jq) and the states with m > doubly- 
occupied sites. 

VI. BIG PICTURE: GENERAL 
TRANSFORMATION 

Thus far, we have shown that while different meth- 
ods of performing perturbation theory preserve the orig- 
inal energy spectrum of the Hubbard model, the effec- 
tive Hamiltonians, Eqs. 10, and Eq. H54|l differ from 
each other and from Eqs. lO, that agree with each 
other. The terms in question all arise at fourth order and 
are all reducible^^ with respect to the zero-double occu- 
pancy sector. That is, they contain hopping processes 
that do not originate from excitation to the doubly oc- 
cupied subspace, TqT_iTiTq, nor terminate once an elec- 
tron is returned to the singly occupied subspace, for ex- 
ample, r_iTiT_iTi. All such processes can be viewed as 



arising from a transformation^- of the eigenstates in the 
low-energy sector. To lay plain how the effective Hamil- 
tonian is unavoidably affected by this transformation, we 
now formulate a general method which makes it possi- 
ble to derive all of the Hamiltonians presented thus far 
within a single computation scheme. Our starting point 
is BW integral equation, Eq. (|14|l . whose solution we 
write symbolically as 

\^)^G{t,E)\^jo). (55) 

The exact expression for the energy-dependent operator 
G{t, E) is obtained by iterating Eq. l(Tl|l . Applying P — 
1—Q (see Sec. I1V|I to the left-hand side of the Schrodinger 
equation, Eq. (|12|l . we obtain 

{E~ToMo)~PT-M)=0, (56) 

which can be recast as a non-linear eigenvalue problem 

S|Vo) = (^0 + PT_iG(t, ^))|^o), (57) 

using Eq. (|55|l . Taylor expansion of G{t,E) results in 
a polynomial in the energy eigenvalue. Through fourth 
order we find that. 



E\yjo) = i^To - T-iTi + T^iToTi - T-iT^Ti - -T\T^ + - ET^iTi - E^T^^Ti + 2ET^iToTij j |^o). (58) 

To eliminate the energy-dependence on the right-hand side of this equation, we substitute Eq. (|57|l for each occurrence 
of E\iPq) until all the energy dependence has disappeared. The result of this procedure is an eigenvalue problem 



(59) 



r 



with a non-hermitian operator H. All the terms in the 
second parenthesis arise explicitly from eliminating the 
energy dependence in Eq. lf5S|) and as a consequence are 
products of the proper BW terms in the first parenthe- 
sis. It is in this sense that such terms are reducible with 
respect to the zero double occupancy sector. The lack of 
hermiticity surfaces because projection does not respect 
the mutual orthogonality of the eigenstates in any of the 
degenerate subspaces. 

Hermiticity can be restored by a suitable global trans- 
formation of the eigenstates within each degeneracy sub- 
space. To proceed, we consider the operator 

Z = l + aT-iTi + bT-iTiTo + cTqT^iTi + dT-iToTi, 

which is explicitly not unitary. As such, it can be used 
to construct an effective Hamiltonian 



by placing appropriate conditions on the coefficients, a, 
6, c, and d so that Ti.cS is hermitian. Through fourth 
order, hermiticity is restored by demanding that 



(61) 



(60) 



The resultant effective Hamiltonian, 
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+ (T_iTi)2 - i(T_iTir2 + T^TiT^i) + -,i2ToT^iTiTo - T^iT^T^ - T^T^^T^) (62) 

I 



contains three reducible terms whose magnitude is set by 
an arbitrary constant 7. These terms are given precisely 
by the "additional" canonical transformation in Eq. (|10|) . 
Because all the terms controlled by the magnitude of 7 
are reducible, they provide no more than a transforma- 
tion of the eigenstates within the degeneracy subspace 
with zero double occupancy. The multiplicity of Hamilto- 
nians we have derived here all arise from different choices 
for 7. For example, within the canonical transformation 
method of Ref. |a, we have 7 = 0. Effective Hamil- 
tonians within a sector with a fixed number of doubly 
occupied sites can only be determined up to an arbitrary 
rotation of the eigenstates within the degeneracy space. 
To understand what happens in the case of BW pertur- 
bation theory, we recall that, although one starts from 
a set of wave functions that are orthogonal in the full 
Hilbert space, the projection into a degeneracy subspace 
is a process that does not respect the orthogonality. To 
get a hermitian Hamiltonian we have to perform a gen- 
eral transformation on the degeneracy subspace in such 
a way that the orthogonality of the projected compo- 
nents is restored. This is what the general non-unitary 
transformation defined by Eq. is doing. Canoni- 

cal transformations, on the other hand, are by definition 
unitary and hence no additional orthogonalization trans- 
formations are necessary^iS 

VII. SOME COMMENTS ON EXPERIMENT 

Several experimental groups^ have pointed out that 
high-temperature superconductors show spectral weight 
rearrangements over the Mott scales^i Rearrangement of 
spectral weight over large energy scales is expected in 
strongly correlated systems simply because many of the 
eigenstates are localized or almost localized. Using the 
Lehman representation, one can easily see that the mo- 
mentum eigenstates probed by photoemission or optical 
spectroscopy, for example, have non- vanishing projection 
on essentially all the true eigenstates of the interacting 
problem. When the Mott gap is closed, this means that 
spectral weight changes will occur over all the energy 
scales when doping or temperature is changed. When 
the Mott gap is opened, this will continue to be the case 
but, nevertheless, the spectral rearrangements over the 
lower Hubbard band will be describable to a high degree 
of accuracy using only the effective low-energy theory, as 
long as one uses the operators that are appropriate for 
the low-energy sector. These operators take into account 
rearrangements in the upper Hubbard band through vir- 



tual states. For example, upon doping by an amount x, 
exact calculations on the Hubbard model show that the 
spectral weight transfered from the upper to the lower 
Hubbard band exceeds 2a;, while 2x was argued to be 
the prediction of the t — J model by some early worki^ 
Even at the level of the t — J model, however, there is a 
correction^iS to the low-energy spectral weight (LESW) 
from the transitions across the Hubbard gap that arise 
from transforming the electron operators. Through order 
(tjUY ^ the LESW agrees well with the exact diagonal- 
ization on small systems. 

VIII. CONCLUSION 

We have studied several methods for performing de- 
generate perturbation theory. We have shown that, to 
fourth order, they lead to low-energy effective theories 
that appear different but, in fact, are all related through a 
unitary transformation in the low-energy subspace, Eqs. 
(|H1), IjlOl) . The necessity of a unitary transformation in 
the low-energy subspace to prove the equivalence of the 
theories does not normally occur in lower-order theories 
and thus is a rather unfamiliar property. The most sys- 
tematic approaches are the two canonical transformation 
methods, the easiest algebraically are the general trans- 
formation (Sec. IVI|I and the resolvent method, while 
the Brillouin-Wigner method, modified a la Raylcigh- 
Schrodinger, becomes rather cumbersome in higher or- 
ders. Nevertheless, the latter method gives some insight 
into the other approaches. In particular, it allows to un- 
derstand the appearance of terms in intermediate states 
that appear, at first glance, to be in the low-energy sub- 
space ("improper terms"). Also, it shows that in the 
low-energy subspace, eigenvectors at order i" have non- 
zero projection on eigenvectors at order only if the 
expansions of the corresponding energies agree to order 
t"~^. The need to orthogonalize and to keep track of the 
order of the energy expansion makes the BW approach in 
practice more delicate to carry out to high order. How- 
ever, the existence of a small parameter makes it com- 
pletely equivalent to the canonical transformation and 
resolvent approaches. 

This work has laid plain how perturbation theory to all 
orders can be formulated unambiguously. Although the 
question of the accuracy of projected schemes in treating 
the full Hubbard model has not been addressed in the 
present work, it is generally assumed that such projected 
schemes should capture all the physics of the original 
Hubbard model to a given order in tjU . Note, however. 
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that this presupposes that in the strong coupHng hmit 
all observables can be expanded in powers of t/U . For 
an infinite system, this may not always be true. Also, as 
is the case with any perturbation theory, such projected 
schemes obviously fail to adequately describe phase tran- 
sitions. For example, at half-filling the interaction tuned 
Mott transitions is beyond the domain of applicability of 
the projected theory. In the case of the doping-induced 
Mott transition, the projected theory cannot address the 
question of the chemical potential, which is of order U . 
Nonetheless, as long as the interaction U is much larger 
than the bandwidth, it is possible to write a unique (up to 
a unitary transformation in the target space) low-energy 
effective theory whose range of applicability is limited 
by the condition that physical observables have a well- 
defined expansion in t/U . We reiterate that consistency 
of such procedures requires that all operators be trans- 
formed as well. Although we presented the results within 
the Hubbard model, all methods are easily generalizable 
to other models. 
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APPENDIX A: MISCELLANEOUS ADDITIONAL 
RESULTS 

The third-order CT generator 5*3, version of Ref. ^ is 
given by 

-Tf r_i + 2Tir_iTi - 2T_irir_i j (Ai) 

-2ToTiTo + 2Tor_iTo. 

The BW third-order expressions for the matrix element 
of an off-diagonal operator (Oi) increasing the number of 
doubly-occupied sites by one is given by 

= - T_,0, + {T_,To - ToT_i)Oi 

+ ^T.iTiT^iOi - ^T\TiOi (A2) 
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